Finite temperature phase diagram of trapped Fermi gases with population imbalance 
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We consider a trapped Fermi gas with population imbalance at finite temperatures and map out 
the detailed phase diagram across a wide Feshbach resonance. We take the Larkin-Ovchinnikov- 
Fuldc-Ferrel (LOFF) state into consideration and minimize the thermodynamical potential to ensure 
stability. Under the local density approximation, we conclude that a stable LOFF state is present 
only on the BCS side of the Feshbach resonance, but not on the BEC side or at unitarity. Further- 
more, even on the BCS side, a LOFF state is restricted at low temperatures and in a small region 
of the trap, which makes a direct observation of LOFF state a challenging task. 

PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 



I. INTRODUCTION 

There has been considerable interest in paired super- 
fluidity of trapped Fermi gases, where the interatomic in- 
teraction can be tuned by varying the external magnetic 
field p], 0, [|[ . Recently, the experimental realization of 
superfluidity in polarized Fermi gases attracts great at- 
tention, where the numbers N-\ and N± of the two atomic 
species undergoing pairing are different [3, [H, [6|. This 
population imbalance is obviously detrimental to super- 
fluidity, since the Cooper pairing requires equal number 
of atoms from both spin components. Therefore, by in- 
creasing the population imbalance from zero, it is ex- 
pected that the BCS pairing state becomes less favorable 
and eventually gives its way to the normal or other exotic 
phases 0, 0, ©, M, d d • 

One of the most fascinating phenomena in unbalanced 
Fermi systems is the Larkin-Ovchinnikov-Fulde-Ferrel 
(LOFF) states, which was first discussed as a "compro- 
mise" candidate exhibiting both pairing and non-zero 
magnetization in the context of superconductors in the 
presence of a magnetic field Q ■ This exotic LOFF state 
is characterized by an order parameter with one or more 
non-zero components at finite momenta q, hence breaks 
translational and rotational symmetry and forms a crys- 
tal of pairing order (i.e., a supersolid). fn the past sev- 
eral decades, the existence of LOFF states was studied in 
various systems fulL including heavy fermions [l4| and 
dense quark matter (lEI |. 

Comparing to the systems mentioned above, ultra-cold 
Fermi gases provide a super clean experimental platform 
with remarkable controllability, so that they can be stud- 
ied with nearly arbitrary interaction strength and pop- 
ulation imbalance. Therefore, after the realization of 
paired superfluidity in resonantly interacting 6 Li atoms 
with population imbalance [J, [jf, the interest of LOFF 
state in these systems has been greatly intensified. Up 
to now, no evidence of LOFF state has been found yet in 
experiments on the polarized Fermi gas. Theoretically, 
some existing studies give significantly different predic- 
tions on the stable region of LOFF state: some works on 
the homogeneous system conclude that at zero tempera- 
ture the LOFF state is confined to a narrow parameter 



region on the BCS side of the resonance, for both cases of 
a narrow and a wide [l7| Feshbach resonance; while 
some other calculations claim a much larger region of 
a stable LOFF state QJ, [H S3], well existing at the 
unitarity point. What is subtle for these calculations is 
to implement a right set of stability conditions, which 
are usually controversial. Considering that there exist 
various possible competing phase configurations for this 
system, one thus needs to carefully distinguish a stable 
phase from some metastable states or unstable saddle 
point solutions. 

In this manuscript, we consider the possibility of LOFF 
state and map out the detailed finite temperature phase 
diagram for fermionic atoms in a trap with population 
imbalance. To ensure stability, we directly minimize 
the thermodynamical potential instead of using the or- 
der parameter equations, as the latter may give un- 
stable or metastable solutions [2l[. The calculation of 
the full landscape of thermodynamical potential sounds 
to us the only method capable to distinguish a local 
metastable configuration from a globally stable phase. 
To deal with the trap, we use the local density approx- 
imation (LDA) which is typically valid unless the trap 
has a stronganisotropy and/or the total atom number 
is small (El l22l [23|. Our main results are as follows, 
(i) We conclude that a stable LOFF state can only be 
present on the BCS side of the Feshbach resonance. On 
the BEC side and at the unitarity, the LOFF state is only 
a metastable state, (ii) Even on the BCS side of the res- 
onance, a globally stable LOFF state is only restricted at 
low temperatures and in a small region of the trap. In a 
wider temperature and spatial region, the LOFF state is 
only a metastable state. Considering experimental limit 
of temperature and resolution, we expect that a direct 
observation of stable LOFF states is challenging for po- 
larized Fermi gases in typical harmonic traps. 

The remainder of this manuscript is organized as fol- 
lows. In Section [Til we discuss our formalism for unbal- 
anced Fermi gases in an isotropic harmonic trap. We 
first derive a Ginzburg-Landau theory for the second 
order normal-LOFF and normal-BCS phase transitions, 
and then derive the mean-field thermodynamical poten- 
tial to study the first order normal-LOFF and LOFF- 
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BCS phase transitions. In the mean-field calculation, for 
simplicity, we focus on the single-plane- wave LOFF state 
(i.e., the FF state). Our main results are presented in 
Section IIII1 where the phase diagrams showing various 
stable or metastable states are illustrated for systems at 
unitarity, and on the BCS and BEC sides of the Feshbach 
resonance. 

II. THE FORMALISM FOR POLARIZED 
FERMI GASES IN A HARMONIC TRAP 

Since the population of closed channel molecules is neg- 
ligible close to a wide Feshbach resonance [13, HH, we 
study a trapped Fermi gas by considering the following 
single-channel Hamiltonian (we use the natural unit such 
that h = ks = 1 throughout this manuscript): 

W = E Ck,<T«k.o- Q k,g- 
k , <7 

+ E y ( k ' k ') a !j/2+k,T a q/2-ka a q/ 2 - k '^ a q/ 2 + k '^'(l) 



k.k'.q 

where and ak,<r are creation and annihilation oper- 
ators for fermions labeled by the spin (hyperfine state) 
indices a —1,1, respectively, and £k,<r — £k ~ /V is the 
fcrmion dispersion ek = k 2 /(2m) shifted by the corre- 
sponding chemical potential. The attractive fermion- 
fermion interaction V^(k, k') can be written in a BCS 
form as V(k, k') = — g, provided that only s-wave con- 
tact interaction is considered. The interaction strength 
g can be connected with scattering parameters through 
the standard renormalization relation 

1 * e£. 



kpa s 



where Ao = mL 3 kF/(4ir) with Fermi momentum fc^ and 
quantization volume L 3 , and a s is the s-wave scattering 
length. In the following discussion, we take the local 
density approximation such that /if = fj,(r) + h, ji^ — 
/i(r) — h, and /i(r) = (j,q — U(r), where U(r) = moj 2 r 2 /2 
represents an external harmonic trap [26] . The chemical 
potential at the trap center and the chemical potential 
imbalance h can be related to the total particle number 
N = iVf + N± and the population imbalance P = (iVj — 
iV;)/(iV T + iV ; ). 

Using the functional integral technique, we introduce 
the standard Hubbard-Stratonovich field A which cou- 
ples to a* a) in order to integrate out the fermions, leading 
to the partition function 



Z = Trl 



= J Z>[At,A] cxp{-S cff [At,A]} 



(3) 



with the effective action 
S c n [At, A] = jf" | £ - Trln (/3G _1 [A])| 



where A = A(q, r) depends on momentum q and imag- 
inary time t, [3 = 1/T is the inverse temperature, and 
G 1 is the inverse fermion propagator 



G- 1 [A(q,r)] = 



-cV+£ T A(q,r) 
At(q,r) -0 r -fj 



(5) 



Up to now we have not incorporated any approxima- 
tion and the effective action Eq. ((4]) is accurate. Next, 
we will discuss in the rest of this section two approxi- 
mation schemes, which can be applied to various situa- 
tions. In section Hi A[ a Ginzburg-Landau (GL) theory 
is presented to study the second order normal-LOFF and 
normal-BCS phase transitions. In section Hi Bi a mean- 
field approach is applied to derive thermodynamical po- 
tential, such that the possibilities of first order normal- 
BCS and LOFF-BCS phase transitions can be analyzed. 



A. Ginzburg— Landau theory and 2nd order phase 
transitions 

In this section, we consider only the possibility of sec- 
ond order phase transitions from normal to superfluid 
phase, including both ordinary BCS and exotic LOFF 
states. Close to the phase transition line, the order pa- 
rameter A(q, t) has small magnitude, hence the effec- 
tive action S e g can be expanded in powers of A. In the 
spirit of Ginzburg-Landau (GL) theory, we are interested 
in static A(q). Therefore, the effective action takes the 
form 



S cff «5>(q)|A(q)| 2 +0(A 4 ), 



(6) 



where a(q) = g^ 1 — xCQjO) with x(q, 0) is the pair sus- 
ceptibility. Notice that as a(q) depends only on the mag- 
nitude of wave vector q, the fourth order terms are cru- 
cial to determine the crystalline structures of the LOFF 
state [l3j]. However, when we getting very close to the 
phase transition line where the effective action is dom- 
inated by the leading quadratic term, it is sufficient to 
consider the coefficient a(q), leading to 



«(q) = E 



i 

2^ 



1 - »F(gkq, + ) - »F(£kq,-) 



2 £kq 



M 

kpa s 



(7) 



where np(x) = (e" x + 1) is the Fermi distribution, 

£kq,± = 6cq±(6e k q-/l), £ k q = J^a (6c+q,t7 + £k,<x )/4, and 

<5ek q = (fk+q — £ k)/2. As concluded in the GL theory, 
the normal state with A(q) = loses its stability as long 
as a(q) becomes negative for one or more q components. 
Therefore, we introduce the order parameter equation, 
which corresponds to the condition that a(q) crosses zero 



(4) 



kpa s 



E 



2fk 2£kq 



(8) 
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FIG. 1: Temperature versus chemical potential phase di- 
agram is illustrated by solving Eq. ((SJl at unitarity, show- 
ing three different regions: (1) normal state is stable or 
metastable, (2) normal state loses its stability at q / 0, 
and (3) normal state loses its stability at q = 0. In this 
figure, the chemical potential difference h = 0.4 with arbi- 
trary energy unit. Inset: The T-r diagram can be obtained 
by incorporating the LDA approximation. Here, we choose 
the energy unit as the chemical potential at the trap center 
fj,o = 1, and h — 0.4 in this scale. The radius r is in unit of 
ro = ^Co/f™ 2 ). 



This order parameter equation can be solved for a given 
li(r) and h to obtain the normal-BCS transition tem- 
perature T c (q = 0), while the normal-LOFF transition 
temperature is determined by the maximal T c (q) for all 
finite q ^ 0. In Fig.[Tl a typical T-fi phase diagram is de- 
picted showing three different regions, where: (1) normal 
state is stable or metastable [a(q) > for all q], (2) nor- 
mal state is unstable due to LOFF instabilities [a(0) > 0, 
but a(q) < for some finite q's], and (3) normal state is 
unstable due to BCS instability [a(0) < 0]. Notice that 
at unitarity, the scattering parameter a s — > oo and does 
not set up a length or energy scale. Therefore, phase 
diagrams for different values of h (as long as h ^ 0) 
are identical when an appropriate energy unit is applied. 
Furthermore, by incorporating the local density approxi- 
mation (LDA) /i(r) = /io — Z7(r), this T-fj, phase diagram 
can be easily transformed into the T-r plane (r is the ra- 
dius of the trap, see inset of Fig. [l}, which can be related 
to experiments. 

It should be emphasized that the phase diagrams ob- 
tained by solving the GL equation © is not a result 
within mean- field [l(| or NSR schemes [27l . [28j . since 
the according approximated number equations are not 
included for self-consistent solving. Therefore, under the 
assumption that the order parameter is small, the GL 
theory discussed above is well controlled and the results 
are reliable. However, the small order parameter restric- 



tion sets a limit of this approach only for analyzing sec- 
ond order phase transitions. The complete phase dia- 
gram, where the possibilities of first order phase transi- 
tions have to be taken into consideration, is beyond the 
scope of GL theory since |A| is not necessarily small. 
Thus, we discuss in the next subsection a mean-field ap- 
proach to derive the thermodynamical potential, which 
offers an approximation technique to study the complete 
phase diagram by direct minimization. 



B. Mean-field theory and 1st order phase 
transitions 

Unlike the GL expansion of effective action Eq. ^ 
discussed above, a mean-field theory is considered here 
by introducing a uniform static saddle point A(q, r) = 
Acj(5(q — Q), which corresponds to the BCS state for 
Q = 0, and to the single-plane- wave LOFF state for finite 
Q. With this assumption, the saddle point action is 



So = P'- 



£{/?(&Q 



Ei 



- In [n F (-E kQi+ )] + In [n F (-E kQt _)} }, (9) 



where E kQ = i/£ kQ + A^, and £ k Q,± = E kQ ± (<5e kQ - 

h) is the quasiparticle energy (+) or the negative of the 
quasihole energy (— ). 

The saddle point conditions SSq/6Aq = and 



SS /5\Q\ 

k F d s 



Q 

lead to the order parameter equations 



E 



_L I - n F (£ , kQ:+ ) - n F (£; kQ ,_) 

2e k 



2£ kQ 



o = EdQi-A^ [i + M£ kQ , + )-n F (£; kQ ,_)] 

k 



E kQ 



[-1 + n F (£; kQ:+ ) + n F (E kQ ,_)] j, 



(10) 



which can be solved together to determine A and |Q|. 
Here, k z = |k|cos0 with 9 is the polar angle. However, 
solving these two order parameter equations is not suffi- 
cient to determine the phase diagram, since the solutions 
may be metastable states or only unstable saddle points. 
Therefore, to ensure stability, we evaluate and directly 
minimize the thermodynamic potential Qq = Sg/0, in- 
stead of imposing various subtle stability criteria [16|, [2(| • 

We show in Fig [5] a typical contour plot of the thermo- 
dynamic potential Qq as functions of A and |Q|, where 
lighter regions denote lower Oo. Notice that the normal 
phase occurs at A = for all values of |Q|. In this 
plot, two local minima are present, corresponding to the 
normal (A = 0) and the BCS (A ^ 0, |Q| =0) states, 
respectively. Combining the results together with those 
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FIG. 2: Contour plots of mean-field thermodynamic potential 
Q,o as a function of A and wave vector |Q|. Two local minima 
(lighter region) are shown to characterize the normal (A = 0) 
and the BCS (A 7^ 0, |Q| =0) states. In this plot, parameters 
are chosen as ft. = 0.4, /1 = 0.5, T — 0.1, and l/(fci?a s ) = 
(at unitarity), with arbitrary energy unit. 

obtained from the GL theory, we can get more infor- 
mation about the finite temperature phase diagrams of 
trapped Fermi gases, which are discussed next. 



III. FINITE TEMPERATURE PHASE 
DIAGRAM OF TRAPPED FERMI GASES WITH 
POPULATION IMBALANCE 

Up to now, we discuss a GL theory which is valid for 
the small order parameter regime but not constrained by 
any approximation schemes, as well as a mean-field ap- 
proach which is approximate but works for wider regions. 
Using these methods and the LDA approximation, we are 
able to analyze the finite temperature phase diagrams of 
polarized Fermi gases in harmonic traps. Next, we first 
consider in Section IIII Al the case of unitarity, where the 
scattering length a s does not set a length or energy scale, 
leading to the presence of universality. In Section IIII B[ 
we discuss the weakly interacting BCS regime, where em- 
phasis is put on the possibility of a stable LOFF state. 
Lastly, the strongly interacting BEC regime is studied in 
Section [mCl 



A. Unitarity 

Following the formalism outlined in the previous sec- 
tion, we first map out the phase diagram for trapped 



fermions at unitarity. The different phases in the trap can 
be identified from the global minimum of the mean-field 
thermodynamic potential at order parameter A and wave 
vector Q. At zero temperature, there are three phases 
which could be possibly present: (i) a BCS superfluid 
state with A ^ at |Q| = 0; (ii) a normal mixed state 
(NM) with A = and two Fermi surfaces (/z j , fjL[ > 0) ; 
and (iii) a normal polarized state (NP) with A = and 
one Fermi surface (/xj > 0, ^1 < 0). The trap bound- 
ary is set when A = and both Fermi surfaces vanish 
< 0). However, at finite temperatures, both the 
trap boundary and the phase boundary between NM and 
NP are blurred. 

We show in Fig. [3Ja) the finite temperature phase di- 
agram for a trapped Fermi gas with population imbal- 
ance. From the trap center to the edge, the BCS, normal 
mixed (NM), and normal polarized (NP) phases are iden- 
tified sequentially. Furthermore, within the BCS phase, 
four regions can be identified due to the existence of nor- 
mal and LOFF metastable states, which can be explicitly 
shown from the thermodynamic potential. At low tem- 
peratures, the thermodynamic potential acquire a double 
well structure near the trap edge. Starting from the trap 
edge, the normal minimum at A = is lower than the 
BCS minimum at (A ^ 0, |Q| — 0), while the relative or- 
der reverses after crossing the normal-BCS phase bound- 
ary [see Fig. |3Jb)]. Continuing towards the trap center, 
the normal minimum becomes unstable due to LOFF in- 
stability, but remains stable for BCS instability. The 
LOFF instability can lead to a metastable LOFF state 
at low temperatures, as shown in Fig. 02c). In this case, 
the global minimum in the landscape still corresponds 
to a BCS state. However, from the normal state to the 
BCS state, one needs to pass a potential barrier through 
a first-order phase transition. Further to the trap cen- 
ter or at higher temperature, the LOFF state loses its 
metastability such that the double well structure disap- 
pears and the BCS state becomes the only minimum in 
the landscape of the thermodynamic potential, as shown 
in Fig.^d). What is special in this case is that to go from 
a normal state to a BCS state (with |Q| =0), one needs 
to follow a path of LOFF instability with pair momentum 
|Q| ^ (there is no BCS instability at |Q| =0 when the 
order parameter A is small). As one moves even further 
towards the trap center, the normal minimum becomes 
unstable due to both LOFF and BCS instabilities, as il- 
lustrated in Fig. OJe), and the system goes to the BCS 
phase directly through the BCS instability at |Q| = 0, 
corresponding to a second order phase transition. 

Considering the universality present at unitarity, the 
qualitative features of phase diagrams are identical for 
arbitrary values of h > 0, as shown in Fig. 2] (the tem- 
perature and phase boundaries get somewhat rescaled). 
Therefore, we can conclude that for a polarized Fermi gas 
in a harmonic trap where the interaction is tuned at res- 
onance, a globally stable LOFF state can not be present 
on the finite temperature phase diagram, although there 
exist a small region of a metastable LOFF state and also 
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FIG. 3: (a) Finite temperature phase diagram of a polarized 
Fermi gas in a trap with h = 0.4 and l/(kFa 3 ) — 0, where 
the chemical potential at the trap center is used as energy 
unit, and ro = \/2/j,o / (mu) 2 ) is used as length unit. The com- 
plete phase diagram is shown in the inset, while the selected 
area (rectangle) is zoomed in to show detailed structures. In 
this plot, a superfluid BCS, a normal mixed (NM), and a 
normal polarized (NP) states can be sequentially identified 
from trap center to trap edge. While the BCS and normal 
regions are separated by a phase transition line (dark solid), 
the trap boundary and the phase boundary between NM and 
NP (dot-dashed) are blurred at finite temperatures. Further- 
more, within the BCS phase, four regions can be identified 
due to the existence of metastable normal and LOFF states 
(gray solid). The contours of thermodynamic potential at rep- 
resentative points of each region are plotted in (b-e) , showing 
corresponding characteristic behaviors. 
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FIG. 4: Finite temperature phase diagrams of polarized Fermi 
gases in a harmonic trap at unitarity l/(kF<is) = 0, where 
the same energy and length units are used as in Fig. |3ja). 
Parameters used in these plots are (a) h — 0.1; (b) h — 0.3; 
and (c) h — 0.5. 
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a region where the normal state becomes unstable due to 
the LOFF instability (but still ends to a BCS state). This 
conclusion is consistent with the findings in Refs. [l(| for 
the zero-temperature case, but does not agree with the 
outcomes in pi [l9l [20j . 



B. BCS regime 

Comparing to the unitarity case, the phase diagrams 
on the BCS side of Feshbach resonances are more com- 
plicated. In Fig. 03a), we show a typical case on the BCS 
regime. The most significant feature is the presence of 
a stable LOFF state at low temperatures, which is char- 
acterized by a global minimum of the thermodynamic 
potential at A(Q) ^ with |Q| ^ 0, as shown in the 
contour plots Fig. [5jf) and (g). 

In addition to the presence of a stable LOFF state, the 
BCS phase also contains four regions which can be identi- 
fied due to the existence of metastable normal and LOFF 
states. By evaluating the thermodynamic potential Oo, 
the four regions can be characterized as Qo\ (i) has two 
minima corresponding to a stable BCS and a metastable 
normal states [see Fig. EJb)]; (ii) has two minima corre- 
sponding to a stable BCS and a metastable LOFF states 
[see Fig.[5]Jc)]; (iii) has only one BCS minimum while the 
normal state is stable for BCS instability, but unstable 
for LOFF instability [see Fig.[S]Jd)]; and (iv) has only one 
BCS minimum and the normal state is unstable for BCS 
instability [see Fig. EJe)]. 

It should be emphasized that unlike the unitarity case, 
universality is not present in the BCS regime, since the 
finite scattering parameter a s sets an additional length 
or energy scale. Therefore, the characteristics of phase 
diagrams for various l/(koa s ) < and h > are not iden- 
tical. However, as shown in a series of phase diagrams in 
Fig. [BJ we can summarize some common features within 
a qualitative level. First, the region where a LOFF state 
is stable is only present in the BCS regime, and disap- 
pears as one moves towards the unitarity [see Fig- E^b)] . 
Second, a stable LOFF state is only present in a small 
region of the trap at low temperatures, and restricted for 
smaller chemical potential difference h. By increasing 
h (or equivalently polarization P), the region for stable 
LOFF state shrinks as shown in Fig.^c). 

Although our calculation results can not be directly 
compared with experiments, where the total particle 
number and polarization are observables instead of /io 
and h, the plots in Fig. [S] outline the general features of 
a finite temperature phase diagram in the BCS regime. 
Therefore, we can conclude that a stable LOFF state may 
be present in harmonically trapped Fermi gases with pop- 
ulation imbalance, but only at low temperatures and in 
a small region of the trap. 




A A 
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FIG. 5: (a) Finite temperature phase diagram for a polarized 
Fermi gas in a harmonic trap, showing the small region where 
a stable LOFF state is present. The complete phase diagram 
is shown in the inset, while the selected area (rectangle) is 
zoomed in to show the detailed structure. The parameters 
used in this plot are h — 0.1 and l/(koa s ) = —0.5 (BCS 
regime), where fco = ^J2mp,o with no is set as the energy unit. 
The radius r is in unit of ro = \/2/j,o / (mu 2 ). Within the 
BCS phase, four regions can be identified due to the stability 
of normal and LOFF states and how they are broken (gray 
solid) . The thermodynamic potential at representative points 
of each region are plotted in (b-f), while (g) shows the detailed 
structure of the selected area (rectangle) in plot (f). 
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FIG. 6: Finite temperature phase diagrams of polarized Fermi 
gases in a harmonic trap in the BCS regime, with the same en- 
ergy and length units as used in Fig.JS^a). In (a) and (b), the 
complete phase diagrams are shown in the insets, while the se- 
lected areas (rectangle) are zoomed in to show detailed struc- 
tures. Notice in these two plots that for fixed chemical poten- 
tial difference h — 0.1, the region with stable LOFF state be- 
comes negligible and the region with metastable LOFF state 
[M-LOFF, also characterized by point c in Fig. [5ja)] shrinks 
as moving towards unitarity. Furthermore, if fixing the value 
of l/(/coo 3 ) = —0.5 as in Fig. [5] the stable LOFF state also 
almost disappears by increasing htoh = 0.3, as shown in (c). 
Other parameters used in these plots are (a) l/(koa s ) = —0.3, 
and (b) l/(fcoa s ) = —0.1, where fco = \f2m[io. 



FIG. 7: Finite temperature phase diagrams of polarized Fermi 
gases in a harmonic trap in the BEC regime. The same energy 
and length units are used as in Fig. [II] Three regions may be 
present in these plots, where the thermodynamical potential 
(1) has only one BCS minimum and the normal state is un- 
stable for BCS and LOFF instabilities; (2) has only one BCS 
minimum and the normal state is unstable only for LOFF 
instability; and (3) has two minima corresponding to a sta- 
ble BCS and a metastable normal states. Parameters used in 
these plots are (a) l/(fcoa s ) = 0.1, h = 0.3; (b) l/(fcoa s ) = 0.3, 
h = 0.3; and (c) l/(fcoa s ) = 0.5, h = 0.3, where ko = y / 27n7io. 
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C. BEC regime 

As concluded in the previous discussion, the LOFF 
state becomes less favorable as one moves from the BCS 
regime to the unitarity. As one further increases the in- 
teraction strength, the same trend is kept in the BEC 
regime, as shown in Fig. [7] Although the finite scatter- 
ing length a s sets an additional energy scale such that 
the phase diagrams are qualitatively different for vari- 
ous l/(fcoa s ) > and h > 0, there are still some gen- 
eral features which can be concluded. First, similar to 
the unitarity stable LOFF state is not present 

on the finite temperature phase diagrams in the BEC 
regime. Second, the structure of the BCS phase be- 
comes simpler and contains less sub-regions. By mov- 
ing from the unitarity towards the BEC limit, the re- 
gion with metastable LOFF state disappears first such 
that at the value of l/(fc a s ) = 0.1 [see Fig. EJa)], only 
three regions are present in the BCS phase. By increasing 
the value of l/(koa s ), the region where the normal state 
becomes unstable only due to LOFF instability [region 
(2) in Fig. 0a)] disappears, such that the BCS phase 
becomes even simpler as shown in Fig. 0b). Further 
towards the BEC side, only a simple BCS superfluid 
state [region (1)] is present at the trap center, leading 
to a phase diagram similar to the unpolarized case [see 
Fig. 0(c)] . Notice that for the parameters used in these 
plots, the breached pair state (BP1) [HHHI characterized 
by gapless fermionic excitations with a Fermi surface is 
not present at zero temperature. At finite temperatures, 
the boundary between BCS and BP1 is blurred and the 
crossover cannot be clearly distinguished. 



[1] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 
[2] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 

Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. 

Lett. 92, 120403 (2004). 
[3] C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. 

Jochim, J. Hecker Denschlag, and R. Grimm, Science 

305, 1128 (2004). 
[4] M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and W. 

Ketterle, Science 311, 492 (2006). 
[5] G. B. Partridge, W. Li, R. I. Kamar, Y. Liao, and R. G. 

Hulet, Science 311, 503 (2006). 
[6] M. W. Zwierlein, C. H. Schunck, A. Schirotzek, and W. 

Ketterle, Nature 442, 54 (2006); Y. Shin, M. W. Zwier- 
lein, C. H. Schunck, A. Schirotzek, and W. Ketterle, 

Phys. Rev. Lett. 97, 030401 (2006); G. B. Partridge, W. 

Li, Y. A. Liao, R. G. Hulet, M. Haque, and H. T. Stoof, 

Phys. Rev. Lett. 97, 190407 (2006). 
[7] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 

20, 762 (1965); P. Fulde and R. A. Ferrel, Phys. Rev. 

135, A550 (1964). 
[8] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 
[9] H. Muther and A. Sedrakian, Phys. Rev. Lett. 88, 252503 

(2002). 



IV. SUMMARY 

We discuss in this manuscript the finite temperature 
phase diagrams of a trapped Fermi gas with popula- 
tion imbalance, focusing on the existence and stability 
of the LOFF state. We first derive a Ginzburg-Landau 
(GL) theory to study the second order normal-BCS and 
normal-LOFF phase transitions, where the order param- 
eter is assumed to be small to ensure validity. Fur- 
thermore, in order to determine the complete phase di- 
agram, we adopt the mean-field approximation and di- 
rectly minimize the resulting thermodynamic potential. 
This method allows us to distinguish the stable, the 
metastable, and the unstable saddle point phases from 
solutions of the order parameter equations. 

Using these methods, we are able to map out the finite- 
temperature phase diagrams over the BCS to BEC re- 
gion. We show that a stable LOFF state exists only on 
the BCS side of the Feshbach resonance, but not at uni- 
tarity or in the BEC regime. Furthermore, we find that 
the LOFF state only exists at low temperatures within 
an appropriate region of population imbalance. With a 
global harmonic trap, even in the most optimal situation 
of all the parameters, the LOFF state is only present in a 
small region of the trap. With the experimental limits on 
temperature and spatial resolution, we expect that it is 
very challenging to make a direct observation of a stable 
LOFF state in typical harmonic traps. 

This work is supported by the MURI, the DARPA, 
the NSF under grant number 0431476, the ARDA under 
ARO contracts, and the A. P. Sloan Foundation. 



[10] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 

(2003) ; 

[11] P. F. Bedaque, H. Caldas, and G. Rupak, Phys. Rev. 

Lett 91, 247002 (2003). 
[12] W. Yi and L.-M. Duan, Phys. Rev. Lett 97, 120401 

(2006). 

[13] R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 

(2004) . 

[14] H. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Han- 
nahs, E. C. Palm, S. W. Tozer, and D. Hall, Nature 425, 
51 (2003). 

[15] M. Alford, J. A. Bowers, and K. Rajagopal, Phys. Rev. 

D 63, 074016 (2001). 
[16] D. E. Sheehy and L. Radzihovsky, Phys. Rev. Lett. 96, 

060401 (2006). 

[17] W. Yi and L.-M. Duan, unpublished, see also the PhD 

thesis of W. Yi. 
[18] K. Machida, T. Mizushima, and M. Ichioka, Phys. Rev. 

Lett. 97, 120407 (2006). 
[19] J. Kinnunen, L. M. Jensen, and P. Torma, Phys. Rev. 

Lett. 96, 110403 (2006). 
[20] Y. He, C.-C. Chien, Q. J. Chen, and K. Levin, Phys. 

Rev. A 75, 021602(R) (2007). 
[21] W. Yi and L.-M. Duan, Phys. Rev. A 73, 031604(R) 



9 



(2006); Phys. Rev. A 74, 013610 (2006). 
[22] T. N. De Silva and E. J. Mueller, Phys. Rev. Lett. 97, 
070402 (2006). 

[23] For a strongly anisotropic trap where the system becomes 
almost one dimensional, the condition for the LOFF state 
may become less stringent, see, e.g., G. Orso, Phys. Rev. 
Lett. 98, 070402 (2007); H. Hu, X.-J. Liu, and P. D. 
Drummond, Phys. Rev. Lett. 98, 070403 (2007). 

[24] Q. J. Chen, J. Stajic, S. N. Tan, and K. Levin, Phys. 



Rep. 412, 1 (2005). 

[25] W. Yi and L.-M. Duan, Phys. Rev. A 73, 063607 (2006). 

[26] We assume here the trap has a spherical symmetry with- 
out loss of generality, as any non-spherical harmonic trap 
can be rescaled to a spherical one under the LDA [21J. 

[27] M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 
Simons, Nature Phys. 3, 124 (2007). 

[28] X.-J. Liu and H. Hu, Europhys. Lett. 74, 574 (2006). 



